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We study both theoretically and experimentally the effects of introducing deliberate disorder in 
a slow-light photonic crystal waveguide on the photon density of states. We first introduce a the¬ 
oretical model that includes both deliberate disorder through statistically moving the hole centres 
in the photonic crystal lattice and intrinsic disorder caused by manufacturing imperfections. We 
demonstrate a disorder-induced mean blueshift and an overall broadening of the photonic density of 
states for various amounts of deliberate disorder. By comparing with measurements from a GaAs 
photonic crystal waveguide, we find good qualitative agreement between theory and experiment 
which highlights the importance of carefully including local field effects for modelling high-index 
contrast perturbations. Our work also demonstrates the importance of using asymmetric dielec¬ 
tric polarizabilities for modelling positive and negative dielectric perturbations when modelling a 
perturbed dielectric interface in photonic crystal platforms. 
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I. INTRODUCTION 

Photonic crystal (PC) cavities and waveguides are 
attractive nanophotonic platforms for controlling and 
studying fundamental light-matter interactions. Aided 
by the presence of a photonic band gap (PEG), which 
arises from the underlying periodic dielectric structure, 
light within a PC cavity or waveguide is strongly con¬ 
fined within a small volume or area. In the case of a PC 
waveguide (PCW), light can be slowed down by orders 
of magnitude compared to a typical slab or ridge waveg¬ 
uide, which increases the local density of photonic states 
(LDOS). The ability to control light-matter interactions 
in PC platforms leads to a host of photonic applications 
and rich optical interactions [1, 2]. For example, PC cavi¬ 
ties have been used for exploring cavity quantum electro¬ 
dynamics (cavity-QED) in both the weak and strong cou¬ 
pling regimes [3, 4], while PCWs have been exploited to 
realize on-chip single photon sources [5-7]. Slow-light in 
PCWs also enhance non-linear processes including pulse 
compression and soliton propagation [8], third-harmonic 
generation [9] and four-wave mixing [10]. In addition, 
PCWs have been integrated in various photonic circuits, 
as optical sensor elements for refractive index measure¬ 
ments in biosensing [11] and chemical fluid detection [12]. 

In practice, PCWs are highly sensitive to manufac¬ 
turing imperfections (intrinsic disorder) which is in¬ 
evitably introduced at the fabrication stage. Disorder- 
induced losses are particularly detrimental in the slow- 
light regime [13], which was predicted theoretically by 
Hughes et al. [14] using a photonic Green function ap¬ 
proach and is now a common finding of various similar 
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theoretical works in the literature [15-18]. With con¬ 
tinued improvements in semiconductor fabrication tech¬ 
niques, and improved theoretical understanding about 
how to mitigate disorder-induced losses [19, 20], various 
groups have experimentally demonstrated PCW designs 
that have reduced disorder-induced losses [21, 22]. 

However disorder in PCs is not necessarily a hindrance 
as was noticed by John in 1987 [23], who proposed disor¬ 
dered PCs for experimentally observing the well known 
phenomena of Anderson localization. Topolancik et al. 

[24] experimentally demonstrated spectral peaks bearing 
signatures of Anderson localization arising from localized 
modes in a deliberately disordered PCW. Patterson et al. 

[25] utilized coupled mode theory to highlight the effect 
of light localization via multiple scattering by examin¬ 
ing the transmission through a PCW in the slow-light 
regime. The strong localized resonances in disordered 
PCWs have also been used to enhance the spontaneous 
emission factor of embedded quantum dots [26], and re¬ 
cently Thyrrestrup et al. [27] have proposed coupling 
quantum dot emitters to a disordered PCW as a promis¬ 
ing platform for conducting QED experiments. Other 
applications of disordered PCWs include enhanced light 
harvesting and random lasing [28]. 

Apart from causing propagation losses and disorder- 
induced localized resonances, disorder also induces 
changes in the eigenfrequencies and eigenmodes of the 
underlying PC. Ramunno and Hughes [29] modelled 
disorder-induced resonance shifts in PC nanocavities and 
predicted a non-trivial disorder-induced mean blueshift 
in the cavity resonance. Patterson and Hughes [30] ex¬ 
tended this formalism to PCWs, and predicted both 
a mean blueshift of resonances and a disorder-induced 
mode edge broadening. To the best of our knowledge, 
this mean blueshift has not been experimentally mea- 
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sured; this is likely because there was no simple ex¬ 
perimental procedure for proving that a mean blueshift 
occurs, especially for an intrinsically disordered PCW. 
Both of the theoretical works mentioned above dealt with 
intrinsic disorder only, which occurs via rapid fluctua¬ 
tions of the air-dielectric interface and highlighted the 
importance of carefully taking into account local fields 
at the interface. Recently, Savona has exploited a guided 
mode expansion technique to compute disorder-induced 
localized modes and the corresponding spectral density, 
which, as expected, showed sharp spectral signatures 
near the mode edge indicative of spatially localized modes 

[31] . 

In this paper, we introduce a model to describe 
disorder-induced resonance shifts and broadening of the 
fundamental mode near the mode-edge which takes into 
account a systematic increase of the disorder parameters 
(i.e., it allows one to model deliberate disorder which can 
be controlled and changed in a systematic way). This 
leads to a predictable trend that causes an increasing 
mean blue shift and a broadening of the photonic density 
of states (DOS). We show how one can extend the theo¬ 
retical models introduced in Refs. [29, 30] to account for 
both intrinsic and deliberate or extrinsic disorder where 
extrinsic disorder is characterized by a deliberate shift of 
the hole centres. We carefully include local field effects 

[32] in our model to compute the first and second-order 
perturbative changes to the eigenfrequencies of the fun¬ 
damental waveguide mode which then allow us to com¬ 
pute the disordered DOS. For disordered PCWs with 
varying extrinsic disorder (see Ref. [33] for details), we 
compute the ensemble averaged DOS via a Monte Carlo 
approach. Experimentally, measurements of vertically 
emitted intensity are taken for GaAs PCW membranes 
with varying amounts of extrinsic disorder. Since the in¬ 
tensity measurements are a direct measure of disordered- 
induced broadening and frequency shift (blueshift) of 
the DOS, we compare our computed DOS with the in¬ 
tensity measurements and the two are found to be in 
good qualitative agreement. The comparison between 
our theory and experimental data demonstrates the im¬ 
portance of including local field effects when comput¬ 
ing disorder-induced changes to the eigenfrequencies and 
eigenmodes of PCWs. While our theory is perturbative, 
the semi-analytical approach is computationally efficient 
and the results offer useful insights in designing disor¬ 
dered PCWs for spontaneous emission enhancements of 
embedded quantum dots. Finally, we also show an exam¬ 
ple of the underlying disorder-induced quasimodes that 
can be obtained on a finite-size PC lattice by comput¬ 
ing the numerically exact Green function using a full 3D 
FDTD approach [6]. 

Our paper is organized as follows. In Section II we 
review our formalism for modelling disorder-induced res¬ 
onance shifts and point out the limitations of some of 
the polarization models commonly used in literature for 
modelling disorder in PGWs. We then introduce our 
extended polarization model for modelling both intrin¬ 


sic and extrinsic disorder and show how it results in a 
non-vanishing first-order frequency shift. In Section III, 
we highlight our approach for computing the disordered 
DOS given disorder-induced resonance shifts. We also 
present a mathematical argument based on photonic 
Green functions that link the disordered DOS to verti¬ 
cally emitted intensity measurements. In Section IV, we 
numerically compute disorder-induced resonance shifts 
and the disordered DOS which we compare with exper¬ 
imental measurements performed on GaAs PGW mem¬ 
branes. While we find qualitatively good behaviour, in 
Section V we discuss some limitations of our perturbative 
model and we also show some numerically exact simula¬ 
tions of finite-size PGWs, which are limited in spatial 
size because of the numerical complexities. In Section 
VI, we summarize the strengths and weaknesses of our 
perturbative semi-analytic approach and discuss our re¬ 
sults in the context of previous reports in the literature. 
We conclude in Section VII. 


II. DISORDER-INDUCED RESONANCE 
SHIFTS AND DISORDER POLARIZATION 
MODELS 


For modelling the effects of disorder on light scatter¬ 
ing in PCWs, we focus our attention on deriving the 
first and second order perturbative changes to the eigen¬ 
frequencies of a dielectric structure. We treat disorder 
as a perturbation and employ perturbation theory tech¬ 
niques adapted to dielectric structures with high index 
contrasts [32, 34]. Denoting the perturbed eigenfrequen¬ 
cies as O' = Wo + where wq is the unperturbed 

eigenfrequency and represents the order per¬ 
turbation, we restrict ourselves to the first and 

second order Aw*^^^ corrections. Since disorder in PCWs 
is statistical in nature, we compute the ensemble average 
over nominally identical disordered PCWs denoted by 
E []. Thus the first-order ensemble averaged correction is 
given by [29] (w dependence is implicit) 


E 




= E[E*(r).P(r)]dr, (1) 

^ dcell 


where E(r) is the unperturbed eigenmode, P(r) is the 
polarization function to characterize the dielectric disor¬ 
der, and the integration is carried out over the primitive 
unit cell of the PC lattice. The fields are normalized 
according to e(r)E* (r) • E (r) dr = I, where £(r) is 
the unperturbed dielectric constant. The ensemble aver¬ 
age of the second-order correction E[Aw(^)], is computed 
similarly [30] 


E 


Aw(2) 


-f //E[E*(r).P(r)E*(r') 


■ P (r')] drdr'. (2) 


To interpret the ensemble averages of the first and 
second order corrections in a statistical sense, we de- 
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note the net frequency shift as Aw (^j , with 

w = Wq + Aw. The net frequency shift Aw is interpreted 
as a random variable whose probability distribution has 
E [Aw(i)] , E [Aw(2)] as its first and second-order mo¬ 
ments, respectively. Hence, given the unperturbed eigen- 
modes, one is left with choosing a suitable polarization 
model to describe the perturbation of the PC lattice. 
Structural disorder in PCWs can be viewed as intro¬ 
ducing additional scattering sites in an otherwise per¬ 
fect PCW lattice. The scattering sites induce dipole mo¬ 
ments resulting in a disorder-induced polarization which 
acts as a source term in the homogeneous Maxwell equa¬ 
tions, thus contributing to scattering of a propagating 
Bloch mode. An alternative picture is that perturba¬ 
tions will disorder the PC band structure and thus the 
DOS, which will result in disorder-induced localization 
modes and scattering in directions that might otherwise 
be forbidden (e.g., a lossless propagating mode will cou¬ 
ple to radiation modes above the light line in the case of 
a PCW slab). 

There are typically two models that have been widely 
used in the photonics community for modelling dielec¬ 
tric perturbations, which denote as weak-index contrast 
Pw(r) and smooth-perturbation Ps(r). The former model 
neglects the problem of field-discontinuities at high- 
index-contrast surfaces, and is well defined for field com¬ 
ponents that are parallel to the interface; in contrast, 
the latter model addresses this field-discontinuity prob¬ 
lem, though it is appropriate for perturbing a surface 
uniformly in a perpendicular direction, e.g., displacing a 
long sidewall in a direction that is perpendicular to the 
wall interface. By way of a simple example, consider a 
simple planar interface between two dielectrics Eq, Es, lo¬ 
cated at r'. When perturbed by a small amplitude Ah, 
the two polarization models are given as: 

Pw(r) = AEA/i(r)E(r)(5(r - r'), (3) 

Ps(r) = AeAh{r) ("E||(r) -h ^^D_L(r)^ 5(r - r'), (4) 

where Ae = Sa — Ss or Ae = eg — Sa depending on the 
direction of the perturbation, i.e. from Eq to Es or vice 
versa, and E||(r), Dj^(r) denote the parallel and per¬ 
pendicular components of the electromagnetic fields rel¬ 
ative to the boundary interface. The weak-index contrast 
model is accurate in systems exhibiting weak-index con¬ 
trast (i.e., |Ae| 1) and is the most popular choice for 
modelling imperfections in dielectric structures such as 
optical waveguides [35]. In high-index contrast systems 
such as PCWs, the quantity Ae|EP is, however, generally 
ill-defined at the interface due to a large step discontinu¬ 
ity in Ej^ [34], hence the smooth-perturbation model is 
likely more appropriate at the interface due to the use of 
continuous field components. The smooth-perturbation 
model is expected to be valid as long as the perturba¬ 
tion is smooth. Both models have been used to compute 
disorder-induced losses in PCWs [14, 16, 30, 36] and have 
yielded a good qualitative understanding of the observed 
disorder-induced loss phenomena. 


If one views the perturbation (smooth or piecewise 
smooth) as introducing scatterers into the system, one 
must take into account their respective polarizabilities 
which in general depend on the direction of the perturba¬ 
tion. The weak-index contrast and smooth-perturbation 
models assign polarizabilities that differ only in sign when 
the direction of perturbation is reversed but remain un¬ 
changed in magnitude. Moreover, the magnitude of the 
polarizability of a scatterer can be drastically different in 
the weak-index approximation as demonstrated by the 
example of a small dielectric sphere in a homogeneous 
background (see Ref. [37]). Therefore, in general for 
piecewise smooth perturbations such as bumps on an in¬ 
terface, it is important to compute polarizabilities that 
correctly take into account the direction of perturbation. 
To address this concern, Johnson et al. [32] introduced 
the bump-perturbation polarization model, denoted by 
Pb(r) to model surface roughness in PCWs as piecewise 
smooth bumps on the interface, where 

Pb(r) = [Eavgai|E||(r)-FE(r) 7 _LD_L(r)] AP(5(r-r'), 

(5) 

where Eavg = and a||, 7 _l denote the polarizabil¬ 

ities (polarizability tensors per unit volume) of the bump 
perturbation and AV is the volume of the disorder bump 
element. This model is valid for arbitrary dielectric con¬ 
trasts and bump shapes, and useful formulas have been 
obtained for rectangular and cylindrical shaped bumps 
[32]. Using the polarizabilities for a cylindrical bump 
shape, this model has been used to model resonance shifts 
caused by intrinsic disorder in PCWs [30] where a mean 
blueshift and broadening of the ideal bandstructure was 
found. As was noted in Ref. [30], resonance shifts in the 
band structure are not predicted by either the weak-index 
contrast or the smooth-perturbation models. 

In this work, we apply the bump-perturbation model 
with cylindrical bump shape polarizabilities, to compute 
disorder-induced resonance shifts in PCWs, and system¬ 
atically investigate what happens with an increase in the 
disorder parameters for shifted holes. We use this model 
to connect to related experiments on deliberately disor¬ 
dered GaAs membranes where embedded quantum dots 
couple to disorder-induced localized modes resulting in 
enhanced spontaneous emission [26, 33]. The PCW we 
consider is a standard WI formed by introducing a line 
defect in a triangular lattice of air holes etched in a semi¬ 
conductor slab Fig. 1(a). The extrinsic disorder pertur¬ 
bation is characterized by a hole centre shift as shown 
in Fig 1(b). The air holes are cylinders so we employ 
cylindrical coordinates (r, 0, z) henceforth. Furthermore, 
the disordered air hole is assumed to have a constant 
cross section throughout the slab thickness. This allows 
us to replace the disorder volume element in Eq. (5) by 
its cross-sectional area AA and the polarizabilities are 
now 2x2 tensors representing polarizability per unit area 
[32]. To first order, the perturbed area A A of the disor¬ 
der element is proportional to |A/i| which quantihes the 
amplitude of the hole centre shift. 




4 


(a) 


•••••••••• 

•••••••••• 

• 99 •• • •• •• 
• 99•••99 99 

9 99 99 9 99 99 
9 99 99 9 99 99 

9999999999 

9999999999 


(b) 



FIG. 1. (Color Online) (a) Schematic of a disordered W1 
with ideal air holes (dark-filled/black circles), disordered 
holes (light/red circles) and the background slab (grey), (b) 
Schematic of a hole centre shift with the centres marked by 
light-filled/orange circles and the direction of perturbation 
given by the solid-black arrow. The magnitude of the per¬ 
turbation is denoted by Ar and the direction is given by (f>. 
Dashed arrows indicate the shifts of the air-slab interface with 
dashed-red/light and dashed-black/dark arrows representing 
positive and negative shifts respectively. 


In light of current experiments studying localization 
modes and resonance shifts as a function of deliberate 
disorder, we extend the disorder model of Ref. 30 to deal 
with both intrinsic and varying extrinsic disorder. While 
the previous model considered rapid radial fluctuations 
of the air-slab interface as the source of intrinsic disor¬ 
der, here we model both intrinsic and extrinsic disorder 
as a net centre shift of the air hole as shown schemat¬ 
ically in Fig. 1(b). Although intrinsic disorder is likely 
best described by rapid radial fluctuations, the choice to 
model intrinsic disorder as a hole centre shift is driven 
by simplicity as one can map rapid radial fluctuations to 
an effective hole centre shift by comparing experimental 
loss data with numerical simulations as demonstrated by 
Garcia et al. [38]; also, the main effect of the disorder 
below is through deliberate disorder. 

Since disorder is stochastic in nature, we denote 
(Ah, (j)) as the random variables quantifying the total 
disorder (extrinsic and intrinsic) in PCWs. The net am¬ 
plitude of the shift |A/i| is constant around the circum¬ 
ference while the sign is determined by the net azimuthal 
direction of the shift, denoted by 4>. The shift of an in¬ 
finitesimally small arc lying on the circular air-slab inter¬ 
face is then given by 


Ah{Ar, (p; 9) 


+Ar, (j) G fl, 
-Ar, 


( 6 ) 


of the intrinsic/extrinsic disorder perturbations, respec¬ 
tively; A/"(ii, cr) denotes a normal distribution with mean 
/i and standard deviation a while U[a,b] denotes a uni¬ 
form distribution on the interval [a, 6]. The net radial 
fluctuation can be broken down into its Cartesian com¬ 
ponents Ax, Ay, which are given below 

Ax = Axi cos{(pi) -I- Are cos{(pe), (7) 

Ay = An sm{(pi) + Ar^ sin((()e). (8) 

The net radial fluctuation is then given as Ar = 
\/ Ax'^ + Ay^ while the net azimuthal direction is sim¬ 
ply (p = tan“^ (^) ■ Comparing to our previous model 

of rapid radial fluctuations [30] , this model lacks the con¬ 
cept of a intra-hole correlation length as all points on the 
hole shift by the same magnitude but in different direc¬ 
tions depending on the angular hole coordinate 6. How¬ 
ever this model is more appropriate for modelling the de¬ 
liberate displacement of the disordered holes performed 
in the experiment. 

To highlight the main difference between the three po¬ 
larization models discussed earlier, let us look at com¬ 
puting the ensemble averaged first-order frequency shift 
E[AwC)] by using Eqs. (3), (4) or (5) in Eq. (1). For 
weak-index contrast and smooth-perturbation models, 
one must compute the expectation of the total disorder 
E[Ah]. If the extrinsic disorder is zero (ue = 0), it is 
trivial to show that E [A/i] = 0 E[Aa;C)] = 0. In 
the case where both intrinsic and extrinsic disorder are 
present, one can still show E[A/i] = 0 as can be ver¬ 
ified via a Monte Carlo simulation. This result is ex¬ 
pected because given any random value for the net ra¬ 
dial displacement Ar, all possible azimuthal directions 
are equally likely; and since we are assigning symmetric 
weights (differing only is sign) to positive and negative 
shifts in these two models, the first order correction van¬ 
ishes. This is in line with previous findings where intrin¬ 
sic disorder was modelled as rapid radial fluctuations [30] . 
However, a non-zero first-order mean frequency shift is 
expected to occur for the bump-perturbation model since 
E[q!||A/i] 7 ^ 0, E[ 7 _LAh] 7 ^ 0. This is because the polariz¬ 
abilities for the shifts that we use in Eq. (5) are asymmet¬ 
ric, i.e. ajj" 7 ^ a[j~, y]*] ^ yj where -b/— denote positive 
and negative shifts respectively. For the second order 
correction E [Aw*^^^], none of the expectation terms van¬ 
ish and therefore the variance of the net mean frequency 
shift is non zero for all three polarizability models. One 
way to test which model is more appropriate is to com¬ 
pare with experiments where the amount of disorder can 
be controlled, and that is precisely what we do below. 


where fl = [9 — ^,9 9 denotes the polar coor¬ 

dinate and Ar quantifies the net radial perturbation. 
A positive bump/shift (-l-Ar) is defined as the air-slab 
boundary shifting into the slab and vice versa for nega¬ 
tive buinps/shifts as illustrated in Fig. 1(b). We denote 
Arj/g ^ |A/"(0, crye)|j 4>i/e ^ W[—TT, tt] as the random vari¬ 
ables for the radial magnitudes and azimuthal directions 


III. DISORDERED DENSITY OF STATES AND 
CONNECTION TO EXPERIMENTS OF 
VERTICAL LIGHT EMISSION 

Since disorder acts to shift and broaden the mode edge, 
a useful quantity for experimental comparison is the DOS 
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p(a;(k)), defined as the number of frequency levels per 
unit volume of k-space. Unlike the concept of bandstruc- 
ture, which is only well-defined in perfectly periodic sys¬ 
tems, the DOS is valid for all structures. It is well known 
that the DOS of an ideal PCW diverges at the mode edge 
since the group velocity vanishes, while for a disorder PC 
structure, the ensemble averaged DOS exhibits a broad¬ 
ened peak around the ideal mode edge where the width 
of the peak is proportional to the amount of disorder 
present in the PC structure [39]. 

To compute the DOS, we first remark that the defini¬ 
tion of DOS bears close resemblance to the mathemat¬ 
ical definition of a probability density function (PDF). 
Hence, just like a histogram generated from a large sam¬ 
ple dataset represents the underlying PDF, the histogram 
generated from a bandstructure represents the DOS. To 
generate a disordered DOS instance, we generate a dis¬ 
ordered bandstructure given by w(fc) = uJo{k) -|- Aaj(fc) 
where the net frequency shift Aui{k) is a random vari¬ 
able assumed to have a normal distribution with mean 
E [Au;^^^(/c)] and variance E [Aa;*'^^(fc)]. The disordered 
bandstructure allows us to calculate an instance of the 
disordered DOS p{uj). One then computes the ensemble 
averaged disordered DOS p{uj) by averaging over many 
such disordered DOS instances. 

Experimentally, the DOS can be obtained by spa¬ 
tially averaging the vertically emitted light intensity mea¬ 
surements in PCWs. To appreciate how the waveguide 
DOS can be measured through vertical emission, con¬ 
sider the waveguide mode Green function without any 
disorder [40]: 

Gwg(r,r';w) = ^[©(a; - 

+ 0(x-x')4*(r')4(r)e-*'=(“-“')], (9) 


onance of the disorder site. Looking at the limit a; —>■ wq, 

Pj(i’d,wo) = Pj(i’d,wo)——i—7-r, (11) 

7 -b Ap,{rd,uJo) 

where we note that the LDOS at the mode edge is no 
longer divergent, and instead pi(rd,a;e) = ^/A (assuming 
Wo = We), which is simply the LDOS from the disordered 
polarizability model. Since this disordered LDOS is now 
connected to light propagation away from the waveguide, 
vertically emitted light will clearly contain signatures of 
the disordered LDOS for the waveguide modes, and thus 
the disordered DOS when spatially integrated. 

An alternative picture of the disordered DOS can 
be obtained by connecting directly to a sum over the 
disordered-induced modes. In PCWs, propagating and 
localized modes couple with radiation modes above the 
light line resulting in vertically emitted intensity. Near 
the mode edge, the DOS increases due to vanishing group 
velocity leading to an increase in the radiation loss rate 
and a broadened peak in the vertically emitted intensity 
spectrum. Other peaks in the spectrum near the mode 
edge indicate the presence of disorder-induced localized 
modes. Given that the vertically emitted intensity is pro¬ 
portional to the radiation loss rate, denoted by 7 , one 
can show that the radiation loss rate is proportional to 
the DOS. Let us assume the disordered quasimodes (or 
quasinormal modes) [42] are known or can be computed, 
denoted by fj(r) where j indexes the quasimodes which 
have complex eigenfrequencies ujj = loj -|- ijj, where the 
quality factor of each resonance is Qj = Wj/ 27 j. Then 
using mode expansion, one obtains the Green function of 
the disordered PGW by [43, 44] 

Gdis(r,r';w) = 7] (12) 


where Q{x — x') is the Heaviside function, k = k{uj) and 
ffc is the ideal Bloch mode. Now consider adding a point 
disorder model, where the disorder causes a polarizability 
with a Lorentzian lineshape (e.g., typical of a disorder- 
induced resonance or an embedded light source such as a 
quantum dot), = Aj (wq — w — 17 ), where wq is the dis¬ 
order induced resonance frequency, 7 is the broadening 
of the resonance, and A is the coupling strength. Using 
a Dyson equation, G = G -b Ga'^G, the Green func¬ 
tion in the presence of the perturbation can be exactly 
obtained through G, where has units of volume (po¬ 
larizability volume) and the Green function has units of 
inverse volume. Defining Pi(rd) = Im[Gii(rd, fd)] where 
Im[] denotes the imaginary component as a measure of 
the (projected) LDOS, we obtain 


/Oi(rd, w) 


I ^ (w-wo)^-b7[7-bApj(rd,w)] 

Pj Td, wj Hpj(rd,a;)]2 ’ 

( 10 ) 


and thus the disordered LDOS now contains signatures 
of the original waveguide LDOS and the underlying res¬ 


and the LDOS of the disordered PCW is 

p(r, w) = —Im [Tr{Gdis(r, r; w)}], (13) 

TTU! 

where Tr[] denotes the trace. From the total LDOS 
one can compute the DOS by integrating over all space. 
Therefore one sees that the radiation loss rate and the 
vertically emitted intensity is inherently linked to the dis¬ 
ordered DOS. Indeed, each one of the underlying quasi¬ 
modes (and every disordered element) has a vertical de¬ 
cay channel associated with vertical decay above the light 
line. 


IV. CALCULATIONS AND MEASUREMENTS 
OF THE DISORDERED-INDUCED RESONANCE 
SHIFTS AND DENSITY OF STATES 


The experimental samples are W1 GaAs membranes 
with a pitch of a = 240 nm and thickness 150 nm with 
an embedded layer of InAs self-assembled quantum dots 
at the centre of the membrane having uniform density 
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of 80 Quantum dots present a very similar re¬ 

fractive index to that of the surrounding membrane ma¬ 
terial (GaAs). In addition, the experiments presented 
in this paper are carried out under high excitation power 
(57 pW m“^) [33], which drives the quantum dots beyond 
saturation, and they become transparent. For these rea¬ 
sons, we consider negligible quantum dot contribution to 
both inelastic and elastic scattering and, thus, we can 
rule out any quantum dot contribution to our model. 
Various samples each measuring 100 pm long are manu¬ 
factured with varying degrees of extrinsic disorder. Ex¬ 
trinsic disorder is introduced via an additional hole cen¬ 
tre displacement characterized by CTe and is varied from 
0.01a = 2.4 nm to 0.05a = 12 nm in 0.01a = 2.4 nm 
steps. The samples are excited and vertically emitted 
intensity is collected as function of wavelength and po¬ 
sition along the waveguide direction /(A, a;) as shown in 
Figs. 3(a,b). The intensity is then spatially integrated 
along the waveguide /(A) = f Idx as shown in Fig. 3(c). 

To connect to these experiments, we model a corre¬ 
sponding W1 PCW (see Fig. 1(a)) with a slab dielectric 
constant suitable for GaAs (e = 12.11), with the follow¬ 
ing parameters: r = 0.295a (hole radius), h = 0.625a 
(slab height). The ideal bandstructure (i.e., with no dis¬ 
order) is plotted in Fig. 2(a), depicting the fundamen¬ 
tal lossless guided mode that spans from 876 nm (there¬ 
after going above the light line) to 930nm (mode edge). 
The intrinsic disorder is approximated as an effective 
hole centre of shift (see Ref. [38]), and is kept fixed at 
(Ji = 0.005a = 1.2 nm. 

The ensemble averaged first and second-order fre¬ 
quency corrections assuming the bump-perturbation 
model for all samples are plotted in Figs. 2(c), (d). The 
expectations in Eqs. (1,2) were computed numerically, 
that is given the statistical parameters for disorder, N 
samples of the set (Ar^, Are, (j)i, (f>e) are drawn from the 
underlying probability distributions which yields N = 
10"^ samples of A/i . The integration is carried out via 
Riemann sums where the step size is chosen to be small 
enough (3 nm) to ensure numerical convergence. 

Interpreting the first and second order frequency cor¬ 
rections as the mean and variance of the net frequency 
shift (see Sec. II), from Fig. 2(c,d), we see that for all 
cases of disorder, the mean frequency shift is a blueshift 
with a standard deviation that increases as the total 
amount of disorder increases. We note that the pre¬ 
diction of a mean blueshift is non-trivial and is solely 
due to the asymmetric polarizabilities present in bump- 
perturbation polarization model. Since one is often con¬ 
cerned with the mode-edge or cutoff of the fundamen¬ 
tal guided mode in PCWs, we find that the mode-edge 
mean shift and variance is roughly of the same order 
as the amount of disorder in the system. For example, 
for extrinsic disorder of 0.02a(4.8nm), the mode edge is 
blueshifted roughly by 4nm with a standard deviation 
of approximately 4nm. From a bandstructure point of 
view, the mean frequency shift and variance act to shift 
and broaden the bandstructure as shown in Fig. 2(b) for 





FIG. 2. (Color online) (a) Photonic bandstructure (in units 
of wave vector versus vacuum wavelength) of the ideal W1 
showing the fundamental (solid) and higher order (dashed) 
guided modes. The light line is shown in black, (b) The 
broadened (blueshifted on average) bandstructure of the un¬ 
perturbed fundamental mode (blue/dark-solid line) for ex¬ 
trinsic disorder of 0.02a(4.8 nm). The three greyscale shades 
indicate the statistical distribution of the perturbed eigenfre- 
quencies that lie within Ter, ±2cr, ±3 (t (dark-grey, medium- 
grey, light-grey) of the unperturbed fundamental mode, where 
a denotes the standard deviation, (c) Ensemble averaged 
first-order eigenvalue corrections for six disordered sam¬ 
ples representing the mean net frequency shift. (d) En¬ 
semble averaged second-order eigenvalue corrections repre¬ 
senting the standard deviation of the mean net frequency 
shift. In both graphs, the intrinsic disorder is kept fixed at 
0.005a(1.2 nm) while the external disorder is varied as follows: 
Oo, (cyan/solid-light-grey), 0.01a(2.4nm) (green/dashed-light- 
grey), 0.02a(4.8nm) (red/solid-medium-grey), 0.03a(7.2nm) 
(magenta/dashed-medium-grey), 0.04a(9.6nm) (blue/solid- 
dark-grey), 0.05a(12nm) (dashed-black). 


the case with 0.02a(4.8nm) of extrinsic disorder. The 
three grayscale shades demonstrate the statistical nature 
of the frequency shift by differentiating between frequen¬ 
cies that lie within wq ± ct, wq i 2(T, wq i 3 (t where a de¬ 
notes the standard deviation and wq is the unperturbed 
frequency. 

The normalized experimental intensity spectra for two 
different amounts of extrinsic disorder along the waveg¬ 
uide is shown in Figs. 3(a,b). Integrating along the 
waveguide direction, the corresponding intensity spectra 
is compared to the ensemble averaged disordered DOS 
for the six samples (considered previously in Fig. 2) in 
Figs. 3(c,d). For now we neglect the contribution of ra¬ 
diation modes to the DOS which scales roughly as 1/A^ 
(Please see Fig. 3(c) for how this might look like). Treat¬ 
ing the DOS as a probability distribution as mentioned 
in Sec. Ill, each DOS instance histogram had a sample 
size of 1000 (number of k-points) and bin resolution of 
0.27nm (200 bins). The ensemble-averaged disordered 
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DOS was calculated from 500 DOS instances. Note, as 
discussed earlier, the DOS at the mode edge formally 
diverges (as the group velocity approaches zero) in the 
absence of disorder but our computed disordered DOS is 
non divergent and shows a pronounced mean blueshift as 
well as broadening caused by the variance of the net fre¬ 
quency shift. This agrees qualitatively well with the ex¬ 
perimental intensity spectra except for the cases of high 
extrinsic disorder cTe = 0.05a(12nm) where the theory 
overestimates the broadening and for ae = 0.02a(4.8nm); 
where the theory predicts a blueshift, but the experimen¬ 
tal intensity spectra is redshifted. The observed redshift 
of the mode-edge is within the computed variance so it is 
either that this discrepancy arises due to the experimen¬ 
tal sample representing only one disorder instance or the 
fabrication method of these particular waveguides where 
a, e.g., proximity effect could introduce an additional un¬ 
known degree of disorder different from the designed one. 


We now discuss some limitations to our perturbative 
approach. For extrinsic disorder values greater than 
or equal to 0.04a(9.6nm), the computed DOS is too 
broad when compared to the measured intensity, see 
Fig. 3(c), (d). This broadening results from the increase 
in standard deviation of mean frequency shift as shown 
in Fig. 2(d). Strictly speaking, our perturbation theory 
computes mode edge resonance shifts and broadening for 
periodically disordered PCWs; that is the primitive unit 
cell is disordered and then repeated indefinitely. This is 
an approximation as in reality the disordered PCW is a 
concatenation of disordered unit cell instances sampled 
from an underlying probability distribution. Moreover 
for the extreme extrinsic disorder case of 0.05a(12nm), 
we can see the signature of new localized mode form¬ 
ing below the mode edge around 945 nm in the intensity 
spectrum (Fig. 3(c)) which our computed DOS cannot 
reproduce since localized modes that form due to cavity¬ 
like defects are naturally not present in a waveguide ex¬ 
hibiting periodic disorder. 


To assess the role of multiple scattering qualitatively, 
we considered incoherent disorder-induced losses in our 
samples, with and without multiple scattering. With the 
mode edge roughly corresponding to a group index of 
rig ~ 50, our computations indicate that for rig > 20, 
we are already in the regime of multiple scattering for 
all amounts of disorder. Therefore, akin to the overesti¬ 
mation of losses without multiple scattering [45] , the pe¬ 
riodic disorder perturbative approach provides an upper 
bound for mode edge broadening and for more realistic 
predictions, a nonperturbative approach is needed that 
takes into account multiple scattering effects. Such an 
approach is very numerically demanding and is beyond 
the scope of this first paper on the topic. However, below 
we show some numerically exact solutions of disorder- 
induced resonances and LDOS for short length PCWs. 
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FIG. 3. (a,b) Experimental normalized intensity spectra ob¬ 
tained by scanning along the waveguide position for two dif¬ 
ferent amounts of extrinsic disorder as indicated in the figure, 
(c) Experimental spatially-integrated intensity for varying de¬ 
grees of extrinsic disorder as labelled in the figure, (d) Calcu¬ 
lated ensemble averaged DOS (normalized) for the fundamen¬ 
tal waveguide mode for the six disordered samples in (c). The 

red-dashed line in the lowest disorder case is given by A 
with A = 0.2 and represents the qualitative contribution of 
radiation modes to the DOS. The amount of intrinsic disorder 
in all samples is cti = 0.005a(1.2nm). 


V. NUMERICALLY COMPUTED DISORDERED 
INSTANCES FROM A FINITE-SIZE PCW 


Having identified the limits of perturbation theory 
above, we now present some brute force calculations of 
the LDOS using full 3D FDTD computations in a dis¬ 
ordered PCW lattice. The numerical complexity is very 
demanding so we are restricted to much smaller waveg¬ 
uide lengths than used in the experiment; also, we can 
only compute a small number of instances which are not 
enough to compute the ensemble average trend shown in 
Fig. 3. This is mainly due to the large memory require¬ 
ments of the simulation volume since it cannot be reduced 
by using symmetric boundary conditions due to symme¬ 
try breaking caused by disorder. Nevertheless, such cal¬ 
culations are useful for getting a physical picture of what 
is happening for a particular instance and section of a 
disordered PCW. 

To show that the DOS varies from instance to in¬ 
stance given the disorder is kept fixed, we calculate the 
projected LDOS p^(r, w) for ten statistically disordered 

























finite-length PCWs, as shown in Fig. 4, by directly com¬ 
puting the numerically exact photonic Green function of 
the PCW (see Ref. 6 for numerical implementation de¬ 
tails) using the 3D FDTD method [46]. The samples 
we simulate are only 7.2 ^m long (30 unit cells). With 
the waveguide cross section in the xy plane, denoting the 
waveguide direction as x and the origin at the centre of 
the waveguide, we compute the LDOS of a y-oriented 
dipole py{0,ui) placed at the anti-node of Cy, which oc¬ 
curs at the origin. The intrinsic/extrinsic disorder values 
are 0.005a(1.2nm), 0.02a(4.8nm), respectively. While 
ten instances which are only 30 unit cells long are not 
enough to conclude the existence of a mean blueshift of 
the mode-edge, the variance in LDOS profiles can par¬ 
tially explain the discrepancy observed in Fig. 3(b) for 
CTe = 0.02a(4.8nm). For completeness. Fig. 4 also shows 
examples of disorder induced localized modes that ap¬ 
pear both above and below the mode edge. These modes 
are formed via multiple scattering in cavity-like defects 
introduced via disorder. 

We highlight that we have found a 2D FDTD method 
to be inadequate for computing the Green function and 
LDOS for the PCW slab. Firstly, the mode edge for 
a 2D PCW with the same structural parameters (apart 
from the slab height) is different (1.2 pm) and secondly, 
a 2D PCW does not possess radiation or leaky modes 
and out-of-plane decay cannot be computed. Hence the 
computed Green function does not accurately capture 
the realistic 3D resonance shifts expected in the LDOS. 
Although a 2D calculation can capture the qualitative 
modal profile of the localized modes, their sensitivity to 
disorder is quite different to 3D quasimodes. Thus one 
requires a 3D FDTD model to compute the LDOS for a 
PCW slab in general. 


VI. DISCUSSION AND CONNECTIONS TO 
PREVIOUS WORKS 

Our theory, though perturbative, provides an intuitive 
and computationally efficient semi-analytical approach to 
producing experimentally relevant results for moderate 
amounts of extrinsic disorder and provides upper bounds 
for high amount of extrinsic disorder. One computation 
which includes computing the ideal Bloch modes, Monte 
Carlo runs for the expectations and Riemann integrals 
for a given amount of disorder takes roughly 3 hours on 
a single-core CPU whereas computing the LDOS of a 
7.2 pm long disordered PCW using 3D FDTD takes ap¬ 
proximately 10 CPU hours for each disorder instance on 
a cluster using 20 multi-core CPU nodes. 

As we have stated before, the bump-perturbation po¬ 
larization model is crucial to our findings. It is not 
the exact shape of the bump that is important (see 
Refs. [29, 32]) but the use of asymmetric polarizabili¬ 
ties that yields a non-zero net mean frequency shift. In 
the context of disorder-induced losses where all three po¬ 
larization models produce similar results, previously we 
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FIG. 4. (Colour online) (Top) Projected LDOS values 
of a j/-oriented dipole centred at the anti-node of e^, com¬ 
puted for ten instances of disorder (blue/dark-solid) with in¬ 
ternal/external disorder values of Ui = 0.005o(1.2 nm), (Je = 
0.02a(4.8 nm) respectively. For reference, the LDOS with 
no extrinsic disorder (red/light-solid) and ideal mode edge 
(dashed-red/light vertical line) are also shown. All LDOS in¬ 
stances are normalized to their own LDOS peak. The length 
of the waveguides was kept fixed at 7.2pm (30 unit cells). 
(Bottom) As highlighted by the black markers (-|-) on the 
Projected LDOS instance, starting from above (right) and 
going below (left) the mode-edge, localized mode intensity of 
is shown. 


have argued that the bump polarization model should be 
best suited for modelling disorder characterized via rapid 
radial fluctuations and the smooth-perturbation model 
should be valid as long as the air-slab interface remains 
nearly circular [30] which is indeed the case considered in 
this work. To resolve this ambiguity, we rely on the com¬ 
parison with experimental findings (see Sec. IV) which 
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indicates that the bump-polarization model (in the ab¬ 
sence of any other models in the literature) is best suited 
for all types of disorder and various disorder-induced phe¬ 
nomena in PCWs. 

Reference 30 highlighted the importance of accounting 
for local field effects in PCWs by computing disorder- 
induced resonance shifts for the three polarization mod¬ 
els mentioned in Sec. II. While the impact of disorder on 
the DOS is qualitatively well known, to our knowledge, 
no one has quantified the expected resonance shifts or 
spectral broadening as a function of disorder and com¬ 
puted the disordered DOS which is found in good qual¬ 
itative agreement with experiments. Previously, spec¬ 
tral broadening of the bandstructure was observed ex¬ 
perimentally by Le Thomas et al. [47] and predicted the¬ 
oretically by Savona [31] whose findings showed increased 
spectral broadening of the bandstructure as the disorder 
increases, but did not predict a blueshift or quantify the 
expected shift or broadening as a function of disorder. 

VII. CONCLUSIONS 

Through theory and experiment, we have shown that 
accurate modelling of local field effects is critical for com¬ 
puting experimentally relevant mean frequency shifts and 
realistic DOS profiles in PCWs. These findings also point 


out the possible limitations of disorder polarization mod¬ 
els that do not include local field effects or include local 
field effects through the use of symmetric polarizabilities. 
For moderate amounts of disorder, our computationally 
efficient semi-analytical perturbative approach yields re¬ 
sults that are in good qualitative agreement with exper¬ 
iments and can be used to compute realistic quantities 
such as Purcell factor enhancements in PCWs with em¬ 
bedded quantum dots. We have also shown examples 
of the numerically exact LDOS for various disordered in¬ 
stances and the underlying disordered-induced resonance 
modes on a small PCW using a rigorous 3D FDTD ap¬ 
proach. Future work will focus on developing a non- 
perturbative approach that takes into account multiple 
scattering to better model high amounts of disorder. 
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